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Abstract 

In this article, we report on the STAR results of open charm hadron production at mid-rapidity 
in p + p and Au+Au collisions at ~\/s NN = 200 GeV and p + p collisions at V* = 500 GeV. 
The measurements cover transverse momentum range from 0.6 to 6 GeV/c for p + p 200 GeV 
collisions, from 1 to 6 GeV/c for p + p 500 GeV collisions and from to 6 Ge V/c for Au+Au 200 
GeV collisions. D° nuclear modification factor and elliptic flow in Au+Au collisions at ysjvw = 
200 GeV are presented. 



1 1. Introduction 

2 The heavy quark production at RHIC is dominated by initial gluon fusion at initial hard 

3 partonic collisions and can be described by perturbative QCD (pQCD) due to their large mass 

4 [TJ. The heavy constituent quark mass is almost exclusively generated through its coupling 

5 to the Higgs field in the electroweak sector, while masses of (u, d, s) quarks are dominated 
e by spontaneous breaking of chiral symmetry (CS) in QCD [2]. This means that charm quarks 
7 remain heavy even if CS is restored, as it likely is in a QGR One expects therefore that charm 
a production total cross section cr^ should scale as a function of number-of-binary-collisions 

9 iVbin- In addition, if charm quarks participate in the collective expansion of the medium, there 

10 must have been enough interactions to easily thermalize light quarks. Hence, charm quark is an 

11 ideal probe to study early dynamics in high-energy nuclear collisions. 

a 2. Analysis Method and Datasets 

13 Invariant yield of charm quark production InvF is calculated as 

InvF, * N * =J- Yi ^ y) **— (1 ) 

Inpj&pj&y Vmg 2np T Ap T Ay BR / &ag . e rec 

14 where N tr i g is the total number of triggered events used for the analysis. Y(pr,y) is the raw 

15 charm hadron signal in each p T bin within a rapidity window Ay = 2. BR is the hadronic 

16 decay branching ratio for the channel of interest. e rec is the reconstruction efficiency including 

17 geometric acceptance, track selection efficiency, PID efficiency, and analysis cut efficiency. /f rag . 
is represents the the ratio of charm quarks hadronized to open charm mesons. And / tlg is the 

19 correction factor to account for the bias between the minimum-bias sample used in this analysis 

20 and the total NSD sample (3). f^g is found to be unity in Au+Au, 0.65 in p + p collisions at ^ = 

21 200 GeV and 0.58 in p + p collisions at V* = 500 GeV. 

22 Y(pj,y) is obtained from fitting the reconstructed invariant mass spectrum (Fig. [TJ) of open 

23 charm mesons through hadronic decays: 
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D°(D°) -> (BR = 3.89%) 

► D°(D^)n ± (BR = 67.7%) 



D 



K-tt+ti* (total BR = 2.63%) 



The identification of daughter particles is done in the STAR experiment [5| at mid-rapidity 
\y\ < 1 at ySjviv = 200 and 500 GeV. The analysis presented herein is done using three datasets; 
the first one collected in year 2009 (iVtrig ~ 105 million 200 GeV p+p collisions), the second one 
collected in 2010 and 2011 (A^g ~ 800 million Au+Au 200 GeV collisions), and the third one 
in 201 1 (N M b ~ 50 million 500 GeV p+p collisions). 

At present, STAR does not have the capability to reconstruct the secondary vertex of D° 
decay; one must calculate the invariant mass of all Kn pairs coming from the vicinity of the 
primary vertex. This results in a large combinatorial background which was reconstructed via 
the mixed-event method (Au+Au dataset), same-charge-sign, and kaon momentum-rotation (p + 
p dataset) and subtracted from invariant mass spectra of all particle pairs |6|. To reconstruct 
D*, one may exploit the softness of D* — > D°n decay; combine low momentum pions with 
D° candidates, i.e. pairs with 1.82 < M(Kn) < 1.9 GeV/c 2 , and plot difference M(Knn) - 
M(Kn) whose resolution is determined by mostly the soft pion high momentum resolution. The 
combinatorial background is reconstructed by side-band (picking Kn pair outside the D° mass 
region) and wrong-sign (picking soft pion with opposite charge) methods. The dominant source 
of systematic uncertainties for both D° and D* analyses is the difference between yields obtained 
from subtractions of combinatorial background from all particle combinations. 




X 2 /ndl 28.4/27 
yield (361.8±26.1)x10 4 
1 .866 ± 0.001 
0.014± 0.001 i 
STAR Preliminary * 

* •«••• 



Au+Au 0.0 < p < 8.0 GeV/c 
_J I I _E I I 



2500 
2000 
1500 
1000 
500 



1.85 1.9 1.95 







/ 

: I I 


* 

— |— right sign 

wrong sign 

side band 



M(Ktc) (GeV/c 2 ) Mass K „(GeV/c 2 ) M(Kroi)-M(Kjt) (GeV/c 2 ) 

Figure 1: Left panel: D° signal in p+p 200 GeV collisions after same-sign background subtraction |3 |. Middle panel: 
D° signal in Au+Au 200 GeV collisions after mixed-event background subtraction. Right panel: D* signal in p + p 200 
GeV collisions with combinatorial background reproduced by wrong-sign and side-band methods (3). 



43 3. Results 

44 3.1. D meson production in p + p collisions 

45 Yields Y(pr,y) are calculated in six pj bins (first two for D°, the next four for D*) in p+p 200 

46 GeV and five pj bins (first for D , the next four for D*) in p+p 500 GeV. The charm cross section 
at mid-rapidity dcr cc /dy was obtained from power-law function []fit to d 2 o- c y(27Tp T dp T dy) = 



n p< 

48 InvF ■ <x N sd, where InvF is obtained from (flj. cr NSU is the total Non-single Diffractive (NSD) 

49 cross section, which is measured at STAR to be 30.0 + 2.4 mb at -\/s — 200 GeV [7|. In the 
so case of -\fs = 500 GeV, there's no STAR measurement yet; cr NSD is extrapolated from 200 
51 GeV measurement with the help of PYTHIA simulation to be 34 mb. The charm production 



1 cFcr" 



I dcr" («-l)(n-2) 

iy < pr >2(„-3)2 



_^PT__\ 

(PTXn-3)) 



2 



10' 



10" 1 



o 
> 

CD 

>I 

T3 
T3 



O-10- 
CM 



10- 7 



Au+Au 200 GeV (D +D°)/2, lyl < 1 
o 0-80% y10 
^» • 0-80%y10+y11 



0-10% [x20] 
T 10-40% [x5] 
A 40-80% [12] 



X 

STAR Preliminary Jp^ * 

♦ p+p D°+D* [x2] * s 

— p+p Levy scaled by (N ) 



E 



-o10 j |e 
^10 4 



10° r 



4 6 8 

p T (GeV/c) 




\li = 500 GeV data 
Ws - 500 GeV power-law fit*' 

FONLL; m =1.5 GeV/c 2 
- FONLL; in=l. 27 GeV/c 2 
] FONLL uncertainty; m =1.5 GeV/c 2 
] FONLL uncertainty; m c =1.27 GeV/c 2 



1 



6 7 8 
p T [GeV/c] 



Figure 2: Left Panel: D° InvY spectra for various centralities, The last four pj bins in p + p collisions are from D + . 
Right Panel: Charm quark production invariant cross section as a function of D meson pr in 500 GeV p+p collisions 
with two FONLL predictions 1 8 1 using normalization and factorization scale equal to charm quark mass m c . 



52 cross section at mid rapidity 



dv 



ly=0 



is 170 ± 45(stat.)!,,(sys.) fib at = 200 GeV and is 



217 ± 86(stat.) ± 73(sys.) fib at V* = 500 GeV. FONLL predictions for p T spectra QD shown in 
Fig.0 

3.2. D production in Au+Au collisions 
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Figure 3: Left Panel: Elliptic flow as a function of pj . Middle panel: D nuclear modification factor Saa as a function 
of pj for most central (blue) and minimum-bias (red) Au+Au collisions with theoretical predictions from two models 
1 1011 111 , Green rectangles around unity represent systematic uncertainties, from left to right, A^in definition uncertainty 
for the most central (2.8%), N\,m definition uncertainty for all Au+Au (7%), and p + p normalization error (8.1%). Right 
panel: The charm production cross section per A'yn as a function of A^in- 

Yields Y(p T ,y) were calculated in eight p T and three centrality bins. dcr^_ N /dy was obtained 
from the integral of 

d 2 ^ l(2np T dp T dy) = InvF ■ cr ine] /N bin (2) 



over pj and is measured to be -j^j ^ = 175 + 13(stat.) + 23(sys.) fib. InvF is obtained from ( 1 

and <x lnel = 42 mb is the total inelastic cross section [9]. To calculate the D° nuclear modification 
factor 7?aa m various centrality bins, we scaled Levy function^jfit to p+p data by iVbini as shown 



2npj dpjdy dv 2miC[nC+mQ{n-2)\ 



ei in the left panel of Fig. [2] and follow the same process for the original power-law function as 

62 discussed in section |3~T| Since enhanced statistics allow more p T bins in Au+Au collisions, we 

63 rely on the extrapolation from the two fits to estimate one source of systematic uncertainty. The 

64 p + p baseline for 7?aa calculation is the arithmetic average of the Levy and the power-law fit 

65 results. The measurement, shown in the middle panel of Fig. [3] reveals strong suppression in 
ee the most central collisions for pj > 2 GeV/c consistent with the prediction of the SUBATECH 
67 group (Gossiaux) model IfTTI and exhibits the maximum of the R& A around p T 1.5 GeV/c. 
ee This agreement with IfTTI might indicate that the maximum is induced by the transverse flow 

69 picked up from the expanding medium through coalescence with light-quarks. 

70 dcr^ N /dy| v= o as a function of TVbin is shown in the right panel of Fig. [5] Within errors, the 

71 results are in agreement and follow the number-of -binary-collisions scaling, which indicates that 

72 charm quark is produced via initial hard scatterings at early stage of the collisions at RHIC. 

73 The FONLL (darker band) and NLO lfl2ll (lighter band) uncertainties are also shown here for 

74 comparison. 

75 In the Left panel of Fig. [3] the measurement of D° elliptic flow V2 is shown. Within large 

76 statistical error bars, D° is consistent with the STAR Non-photonic electrons V2 indicating non 

77 zero elliptic flow of D° mesons in Au+Au collisions at -v/^nn = 200 GeV. 

78 4. CONCLUSIONS 

79 New open charm hadrons (D°,D* + ) measurements in p + p and Au+Au minimum bias 
so collisions at y«NN = 200 GeV from STAR shows the 7V co n scaling of the charm quark pro- 
si duction cross section at mid rapidity. dcr^_ N /dy| v= o = 170 ± 45(stat.)^](sys.) pb in p + p, 

82 175 ± 13(stat.) ± 23(sys.) pb in Au+Au collisions at 200 GeV and 217 + 86(stat.) + 73(sys.) pb 

83 in p + p collisions at 500 GeV. 

84 The new D° nuclear modification factor /?aa measurement reveals strong suppression in the 

85 most central collisions for pj > 2 GeV/c consistent with the prediction of the model IfTTI and 
se exhibits the maximum of the Raa around pj ^1.5 GeV/c. 

87 In the near future the STAR Heavy Flavor Tracker |[T3l will provide the necessary resolution 

se to reconstruct secondary vertices of charm mesons, which will increase the precision of charm 

89 measurements. 
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